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Abstract. The present paper is devoted to the study of the well-posedness issue for the density- 
dependent Euler equations in the whole space. We establish local-in-time results for the Cauchy 
problem pertaining to data in the Besov spaces embedded in the set of Lipschitz functions, 
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including the borderline case B ^ (R ). A continuation criterion in the spirit of the celebrated 
one by Beale-Kato-Majda in [3] for the classical Euler equations, is also proved. 

In contrast with the previous work dedicated to this system in the whole space, our approach 
is not restricted to the L 2 framework or to small perturbations of a constant density state: we 



Qh ' just need the density to be bounded away from zero. The key to that improvement is a new a 

priori estimate in Besov spaces for an elliptic equation with nonconstant coefficients. 

The evolution of the density p = p(t, x) £ M. + and of the velocity field u = u(t, x) € M. N of a 
nonhomogeneous incompressible fluid satisfies the following density- dependent Euler equations: 

dtp + u ■ Vp = 0, 

p(d t u + u ■ v«) + vn = pf, 

div u = 0. 



>: a) 



Above, / stands for a given body force and the gradient of the pressure VII is the Lagrangian 
\ multiplier associated to the divergence free constraint over the velocity. We assume the space 

■ variable x to belong to the whole W N with N >2. 



A plethoric number of recent mathematical works have been devoted to the study of the 
classical incompressible Euler equations 

f d t u + u ■ Vu + vn = /, 



X ■ (2) 

5_i ■ I divu = 0. 



which may be seen as a special case of ([T]) (just take p = 1). 

In contrast, not so many works have been devoted to the study of (pQ) in the nonconstant 
density case. In the situation where the equations are considered in a suitably smooth bounded 
domain of M? or 1R 3 , the local well-posedness issue has been investigated by H. Beirao da Veiga 
and A. Valli in [H [51 [6] for data with high enough Holder regularity. The case of data with 
W 2,p regularity has been studied by A. Valli and W. Zajaczkowski in [19] and by S. Itoh and A. 
Tani in [TJj. The whole space case 1R 3 has been addressed by S. Itoh in [13J. There, the local 
existence for initial data (pq,uq) such that po is bounded, bounded away from and such that 
Vpo £ H 2 , in uq is in H 3 has been obtained. In the recent paper [TT|, we have generalized 
|13j's result to any dimension N > 2 and any Sobolev space H s with s > 1 + N/2. Data in the 
-+i 

limit Besov space -B 2 2 i are also considered. 

Let us also mention that, according to the work by J. Marsden in [16], the finite energy 
solutions to ([1]) may be interpreted in terms of the action of geodesies. This latter study is 
motivated by the fact that, as in the homogeneous situation, System ([1]) has a conserved energy, 
namely 

(3) \\(VPu)(t)\\h = WVPo^oWh +2 f [ (pf-u)(r,x)drdx. 
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Motivated by the fact that, in real life, a fluid is hardly homogeneous, we here want to study 
whether the classical results for homogeneous fluids remain true in the nonhomogeneous frame- 
work. More precisely, we aim at investigating the existence and uniqueness issue in the whole 
space and in the LP framework for densities which may be large perturbations of a constant 
function: we only require the density to be bounded and bounded away from zero and to have 
enough regularity. We shall also establish blow-up criteria in the spirit of the celebrated one by 
Beale-Kato-Majda criterion for ([2]) (see [3]). 

The functional framework that we shall adopt - Besov spaces embedded in the set C 0,1 of 
bounded globally Lipschitz functions - is motivated by the fact that the density and velocity 
equations of ([T|) are transport equations by the velocity field. Hence no gain of smoothness may 
be expected during the evolution and conserving the initial regularity requires the velocity field 
to be at least locally Lipschitz with respect to the space variable. In fact, as regards the velocity 
field, the spaces that we shall use are exactly those that are suitable for ([2]) . We thus believe 
our results to be optimal in terms of regularity. 

Now, compared to the classical Euler equations, handling the gradient of the pressure is much 
more involved. To eliminate the pressure, the natural strategy consists in solving the elliptic 
equation 

div (aVIl) = div F with F := div (/ — u ■ Vit) and a:=l/p. 
If a is a small perturbation of a constant function a then the above equation may be rewritten 

aAU = div ((a - a) VII) + div F. 

Now, if 1 < p < oo then the standard LP elliptic estimate may be used for absorbing the first 
term in the right-hand side. Hence we expect to get the same well-posedness results as for ([2]) in 
this situation. As a matter of fact, this strategy has been successfully implemented by Y. Zhou 
in [22]. 

In the general case of large perturbations of a constant density state, solving the above 
equation in the l w framework for F G IP may be a problem (unless p = 2 of course). In 
fact, to our knowledge, even if a is smooth, bounded and bounded away from zero, there is no 
solution operator TC : F — > Vn such that 

||vn|| LP < c\\f\\ lp 

unless p is "close" to 2 (see the work by N. Meyers in [IT]). However that closeness is strongly 
related to whether a itself is close to a constant hence no result for all p may be obtained by 
taking advantage of Meyers' result. 

In the present work, we shall overcome this difficulty by requiring the data to satisfy a finite 
energy condition so as to ensure that F is in I? . Indeed, this will enable us to use the classical 
L 2 estimate and from that, it turns out to be possible to get estimates in high order Besov 
spaces Bp r . 

This observation is the conducting thread leading to the first two well-posedness results stated 
in the next section. The rest of the paper unfolds as follows. In section [2j we introduce the 
Littlewood-Paley decomposition and recall the definition of the nonhomogeneous Besov spaces 
Bp r . Then, we define the paraproduct and remainder operators and state a few classical results 
in Fourier analysis. Section [3] is devoted to the proof of existence results and a priori estimates 
for an elliptic equation with nonconstant coefficients in the Besov space framework. To our 
knowledge, most of the results that are presented therein are new. Sections [H [5] and [6] are 
dedicated to the proof of our main existence and continuation results. Some technical lemmas 
have been postponed in the appendix. 

Notation. Throughout the paper, C stands for a harmless "constant" whose exact meaning 
depends on the context. 

For all Banach space X and interval / of R, we denote by C(I;X) (resp. Cb(I;X)) the set 
of continuous (resp. continuous bounded) functions on / with values in X. If X has predual 

v* j-i J j-_ i n It. v\ _ _n ]_J r, j.: £ . t v i_ 
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for any <fi G X* , the function t i— ► (f(t),4>)xxx* is continuous over /. For p G [l,oo], the 
notation L P (I]X) stands for the set of measurable functions on / with values in X such that 
t I— ► belongs to L P (I). We denote by L P oc (I) the set of those functions defined on / 

and valued in X which, restricted to any compact subset J of I, are in L P {J). 
Finally, for any real valued function a over R , we denote 

a* := inf a(x) and a* := sup a(x). 
1. Main results 

As explained in the introduction, we hardly expect to get any well-posedness result if the 
initial velocity is not in C ' 1 . It is well-known (see e.g. [2], Chap. 2) that the nonhomogeneous 
Besov space Bp r is continuously embedded in C ' 1 is and only if the triplet (s,p, r) Elx [1, oo] 2 
satisfies the following condition: 

(C) s > 1 + N/p or s > 1 + N/p and r = I. 

This motivates the following statement concerning the existence of smooth solutions with finite 
energy: 

Theorem 1. Let (s,p,r) satisfy Condition (C) with 1 < p < oo. Let uq be a divergence- 
free vector-field with coefficients in L 2 n B^ r . Suppose that the body force f has coefficients in 
To, To]; Bp r ) n C([— Tq, To]; L 2 ) for some To > 0. Assume that po is positive, bounded and 
bounded away from zero and that Vpo G Bp~ . If p < 2, suppose in addition that (pa — p) G T p * 
with p* := 2p/(2 — p) for some positive real number p. 

There exists a time T £ (0, To] such that System ([1]) supplemented with initial data (po, uo) 
has a unique local solution (p, u, VII) on [— T, T] x with: 

• p ±l e C b ([-T,T] x R N ), Dp G ^([-^T];^- 1 ) fand (p - p) G C([— T, T];L P ") if 
p<2), 

• n G C 1 ( [-T, T] ; L 2 ) n C w ( [-T, T] ; B» r ) , 
. Vn G C([-T,T]; L 2 ) n L\[-T, T\;B' jT ). 

Besides, the energy equality ([3]) is satisfied for all t G [— T, T], and time continuity holds with 
respect to the strong topology, if r < oo. 

A few comments are in order: 

• For the classical incompressible Euler equations ([21), the above result statement (without 
the L 2 assumption) belongs to the mathematical folklore. It has been established in e.g. 
|21j in the case 1 < p < oo and in e.g. [2], Chap. 7 in the case 1 < p < oo. 

-+i 

• Note that the above statement covers the borderline case B p p 1 for Condition (C) 
without any smallness assumption. Thus, up to the lower order L 2 assumption which 
is needed to control the low frequencies of the pressure, it extends the result [22] by Y. 
Zhou mentioned in the introduction. 

• If one makes the stronger assumption that (po — p) G Bp r for some positive constant p 
then we get in addition (p — ~p) G C([—T,T];B^ r ) (or C w ([—T,T];B^ r ) if r = oo). 

• If 1 < p < 2 then no G Bp r implies that no G L 2 . Furthermore, in dimension N > 3, 
the assumption that (po — p) G LP may be omitted if p > N/ (N — 1). Therefore, except 
if N = 2 and p < 2 or if N > 3 and p < N/(N — 1), the density need not to tend to 
some constant at infinity. 

• In contrast with the homogeneous case, in dimension N = 2, the global well-posedness 
issue for ([T]) with nonconstant density is an open (and challenging) problem. Indeed, the 
vorticity u> := d\u 2 — c^n 1 satisfies 

d t u + u ■ Vu + di(±)d 2 LI - d 2 (j;)d 1 U = 0, 
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Under the assumptions of Theorem[TJ the solutions to ([!]) satisfy the following Beale-Kato-Majda 
type continuation criterion. For simplicity, we state the result for positive times only. 

Theorem 2. Under the hypotheses of Theorem {J\ consider a solution (p, u, VII) to ([I]) on 
[0, T) x with the properties described in Theorem [7J If in addition 

(4) / (||Vit||z,°o + ||VII|| B s-i) dt < oo 

Jo p ' r 

then (p, u, VII) may be continued beyond T into a solution of ([T]) with the same regularity. 

Moreover, in the case s > 1 + N/p, the term Vm may be replaced by curl-u in @. 

Remark 1. The above statement has two important consequences: 

• First, as Condition (C) implies that Bp~ l is embedded in L°°, one can show by means 
of an easy bootstrap argument that for data in Bp r , the lifespan of a solution in Bp r 

is the same as the lifespan in B p 1 (which is the larger space in this scale satisfying 
Condition (C) ). 

• Second, by combining the previous remark with an induction argument, we see that if we 
start with smooth data (ao,«o) such that the derivatives at any order of Vao an d uo 
are in in LP then we get a local-in-time smooth solution with the same properties. In 

addition, as above, the lifespan for that smooth solution is only determined by the B p p 1 
regularity. This generalizes prior results in the Holder spaces framework in the case of 
a bounded domain obtained in [6]. 

In the two-dimensional case, the assumption that uq G I? is somewhat restrictive since if, say, 
the initial vorticity is in the Schwartz class then uo G I? implies that the vorticity has average 
over M. N . This motivates the following statement which allows for any suitably smooth initial 
vector-field with compactly supported vorticity. 

Theorem 3. Let Tq be in ]0,oo[ and let (s,p,r) satisfy Condition (C) with 2 < p < 4. Let uq 
be a divergence- free vector-field with coefficients in Bp r . Assume that pq is positive, bounded 
and bounded away from zero, and that X7p G BL* ~ . Finally, suppose that the body force f has 
coefficients in L 1 ([To, To]; Bp r ) and that the potential part Qf of f is in C([—Tq,Tq];L 2 ). 

There exists a time T G (0, To] such that System ([I]) supplemented with initial data (po,uq) 
has a unique local solution (p, u, VII) on [— T, T] x with: 

• p ±x eC b ([-T,T] xR N ), Dp G C W ([-T,T\; B^ 1 ), 

• ueC w ([-T,T];B^ r ), 

. Vn G C([-T,T];L 2 ) n L\[-T,T\;Bl r ). 
Besides, time continuity holds with respect to the strong topology, if r < oo, and the continuation 
criteria stated in Theorem also hold under the above assumptions. 

As regards the well-posedness theory, the study of the limit case p = oo is of interest for 
different reasons. First, the Besov space B^ 1 is the largest one for which Condition (C) holds. 
Second, the usual Holder spaces belong to the family B^ r (take r = oo) and are suitable for 
the study of the propagation of tangential regularity in (fTJ) , and of vortex patches than we plan 
to do in future works. 

Theorem 4. Assume that u G B^ r n LP, f G L 1 ([-T ,To];-B^ o r . n L p ) and p G B^ r for 
some p G (l,oo) and some s > 1 (or s > 1 if r = 1). There exists a constant a > depending 
only on s and N such that if, for some positive real number ~p we have 

(5) ||po -p\\b^ t < oip, 

then there exists some T > such that System §Q) has a unique solution (p, u, VII) with 
. peC{[-T,T];B^ r ) (orC w ([-T,T];B^ r ) ifr = oo), 

_ „. r- nl\ T 1 Ti\. T V r\ IDS \ / „. ^ P I \ T TH . T V \ r*> H l\ T Tl. OS \ Xt ~„ ) 
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. vnei^-T.Ti^y. 

Note that our result holds for small perturbations of a constant density state only. The 
reason why is that, in contrast with the previous statements, here bounding the pressure relies 
on estimates for the ordinary Laplace operator A. In other words, the heterogeneity (a — a) VII 
is treated as a small perturbation term. We expect this smallness assumption to be just a 
technical artifact. However, removing it goes beyond the scope of this paper. 

2. Tools 

Our results mostly rely on the use of a nonhomogeneous dyadic partition of unity with respect 
to the Fourier variable, the so-called Littlewood-Paley decomposition. More precisely, fix a 
smooth radial function % supported in (say) the ball B(0, 3), equals to 1 in a neighborhood of 
B(Q, |) and such that r 1— ► x(re r ) is nondecreasing over M + , and set <p(£) = ~~ 
The dyadic blocks (A g ) gg ^ are defined b}0 

A q := if q < -2, A_i := x(D) and A q := f(2~ q D) if q > 0. 
We also introduce the following low frequency cut-off: 

S q u:= X (2- q D)= £ A p for q > 0. 

p>q-l 

The following classical properties will be used freely throughout in the paper: 

• for any u £ S', the equality u = J2qA q u makes sense in S' ; 

• for all u and v in S', the sequence {S q -\u A q v) q ^ is spectrally supported in dyadic 
annuli. Indeed, as Suppx C -8(0, f) and Supp<£ C {£ 6 M n / § < |f | < §}, we have 

Supp {F{S q -iu A q v)) C e R N I i • 2<? < |e| < f ■ 2«}. 
One can now define what a Besov space S| r is: 

Definition 1. Let u be a tempered distribution, s a real number, and 1 <p,r< 00. We set 

1 

\\u\\ B s r := (^22 rqs \\A q u\\ r LP ^ if r < 00 and ||«||b| iOC := sup 2 9S || A q u\\ LP . 

We then define the space B^ r as the subset of distributions u £ S' such that \\u\\B^ r is finite. 

The Besov spaces have many interesting properties which will be recalled throughout the 
paper whenever they are needed. For the time being, let us just recall that if Condition (C) 
holds true then Bp r is an algebra continuously embedded in the set C 0,1 of bounded Lipschitz 
functions (see e.g. [2], Chap. 2), and that the gradient operator maps B p r in B^~ . The 
following result will be also needed: 

Proposition 1. Let F be a smooth homogeneous function of degree on M. N \ {0}. Then for 
all p £ (l,oo), Operator F(D) is a self-map on LP. In addition, if r G [1, 00] and s£l then 
F{D) is a self-map on Bp r . 

Proof. The continuity on LP stems from the Hormander-Mihlin theorem (see e.g. |12j). The 
rest of the proposition follows from the fact that if u £ B s pr then one may write, owing to 
F(2- q £) = F(f ) for all q > and f ^ 0, 

F(D)u = F(D)A_ lU + ^2(F(p)(2- q D)A q u 

q>0 

where <p is a smooth function with compact support away from the origin and value 1 on the 
support of tp. Note that J-~ 1 (F(p) is in L . Therefore, the standard convolution inequality 
implies that 

\\(Fip)(2- q D)A qU \\ LP <C\\A q u\\ LP 
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while the L p continuity result implies that 

\\F(D)A- X v,\\lp < ||A_i«||£p. 
Putting these two results together entails that F(D) maps Bp r in itself. □ 

Remark 2. Both the Leray projector V over divergence free vector-fields and Q := Id — V 
satisfy the assumptions of the above proposition. Indeed, in Fourier variables, we have for all 
vector-field u with coefficients in S'( 



2O = -j|2^(0- 

The following lemma (referred in what follows as Bernstein's inequalities) describe the way 
derivatives act on spectrally localized functions. 

Lemma 1. Let < r < R. A constant C exists so that, for any nonnegative integer k, any 
couple (p, q) in [1, oo] 2 with q > p > 1 and any function u of L p , we have for all A > 0, 

SuppuC B(0,XR) =► \\D k u\\ Lq < C k+1 X k+N( p'^\\u\\ L p; 
SuppS C {£ 6 R N /rX < |f| < RX} => C- k - 1 X k \\u\\ LP < \\D k u\\ LP < C k+1 X k \\u\\ LP . 

The first Bernstein inequality entails the following embedding result: 

Proposition 2. The space Bt^ r is embedded in the space Bt^ r whenever 

1 < Pi < P2 < 00 an d S2 < si — N/pi + N/p2- 

Remark 3. Recall that for all s£K, the Besov space £>f 2 coincides with the nonhomogeneous 
Sobolev space H s . Furthermore if, for k £ N, we denote by W k ' p the set of L p functions with 
derivatives up to order k in L p then we have the following chain of continuous embedding: 

B k , W k ' p ^ B k 

p,l p,<x 

Let us now recall a few nonlinear estimates in Besov spaces. Formally, any product of two 
tempered distributions u and v, may be decomposed into 

(6) uv = T u v + T v u + R(u, v) 

with 

T u v:=^2S q -iuA q v, T v u :=^2S q -ivA q u and R(u,v) :=^2 A q uA q ,v. 

g g q \q'-g\<i 

The above operator T is called "paraproduct" whereas R is called "remainder" . The decompo- 
sition ([6]) has been introduced by J.-M. Bony in [7]. We shall sometimes use the notation 

T' u v := T u v + R(u,v). 

The paraproduct and remainder operators have many nice continuity properties. The following 
ones will be of constant use in this paper (see the proof in e.g. [2], Chap. 2): 

Proposition 3. For any (s,p,r) £ R x [l,oo] 2 and t < 0, the paraproduct operator T maps 
L°° x B^ r in Bp r , and -B^oo x ^p,r ^ n ^pjr*- Moreover, the following estimates hold: 

\\T u v\\ B s r < C\\u\\ L °°\\Vv\\ B s-i and ||T u T;|| fl .+* < C\\u\\ Bio J\Vv\\ B s-i. 

For any (si,pi,r\) and (s2, P2, ^2) mKx[l,cof such that s\ + S2 > 0, 1/p := 1/pi + I/P2 < 1 
and 1/r := 1/n + I/V2 < 1 the remainder operator R maps B p n iri x BL^ r2 in B^jf S2 . 

Combining the above proposition with Bony's decomposition ([6]), we easily get the following 
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Corollary 1. Let a be a bounded function such that Va G B^ 1 for some s > and (p,r) G 
[1, oo] 2 . T/ten /or any 6 G Bp r n L°° we /ia?;e ao G -Bp r n L°° and there exists a constant C 
depending only on N, p and s such that 

\\ a b\\B* r < \\a\\L°°\\b\\B° T + IHU°°||-Da||Rs-i- 

Pi' Pi I ^PjV 

The following result pertaining to the composition of functions in Besov spaces will be needed 
for estimating the reciprocal of the density. 

Proposition 4. Let I be a bounded interval of R and F : L R a smooth function. Then for 
all compact subset J C L, s > and (p,r) G [l,oo] 2 there exists a constant C such that for all 
a G Bp r with values in J, we have F(a) G B s p r and 

\\F(a)\\ B s r <C\\a\\ B s r . 

Our results concerning Equations (pQ) rely strongly on a priori estimates in Besov spaces for 
the transport equation 

f d t f + v • V/ = g, 
I f\t=o = fo- 



(T) 



We shall often use the following result, the proof of which may be found in e.g. [2], Chap. 3, or 
in the appendix of [8]. 

Proposition 5. Let 1 < p, r < oo and a > 0. Let fo G Bp r , g G L l {[Q,T]; B^ r ) and v be a 
time dependent vector-field in C&([0, T] x R w ) such that for some p\ > p, we have 

N_ 

Vv G L 1 {[{),T];B^ 00 r\L 00 ) if a < 1 + 

Vv G LiQCT];^- 1 ) if <r>l+£, or a = 1+g and r = 1. 

T/ien Equation (T) /ias a unique solution f in 

• the space C([0,T]; Bp r ) if r < oo, 

• fte space (fV <CT C([0, T]; f| ^([0, T]; B^) if r = oo. 



Moreover, for all t G [0, T], we /iai>e 

(7) e-^W||/(t)|| BpV < ||/ || BpV + f e- cv ^\\g{t')\\ BZr dt' 



f ||V«(t)[| * <r<l + £, 

with V'(t) := { B# i00 ni/» 

{ \\Vv(t)\\ K - } if a > 1 + or a = 1 + f ana 1 r = 1. 

If f = v then, for all a > 0, Estimate Q ZioWs wrai/i V'(i) := ||Vu(£)||z,°°. 

3. Elliptic estimates 

In this section, we want to prove high regularity estimates in Besov spaces for the following 
elliptic equation 

(8) - div (oVII) = div F in R^ 
where a = a(x) is a given suitably smooth bounded function satisfying 

(9) a* := inf a(x) > 0. 

Let us recall that in the case a = 1 the following result is available: 

Proposition 6. If a = 1 and p G (1, oo) then there exists a solution map F VLT continuous 
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Proof. We set VII = V(— A) _1 divF. Obviously the pseudo-differential operator V(— A) _1 div 
satisfies the conditions of Proposition [TJ Hence F i— > VII is a continuous self- map on LP . □ 

We now turn to the study of (|8]) for nonconstant coefficients. For the convenience of the 
reader let us first establish the following classical result pertaining to the L 2 case. 

Lemma 2. For all vector-field F with coefficients in I? , there exists a tempered distribution IT, 
unique up to constant functions, such that VII G L 2 and Equation (J8]) is satisfied. In addition, 
we have 

(10) a*||VII|| L 2 < ||F|| L2 . 

Proof. The existence part of the statement is a consequence of the Lax-Milgram theorem. Indeed, 
for A > 0, consider the following bilinear map: 

b\(u,v) = {aVu\Vv)i2 + \(u\v)i i 2 for u and v in H 1 (1S. N ). 

Obviously b\ is continuous and coercive, hence, given F G (L 2 (W N )) N , there exists a unique 
II A G H 1 ^") so that 

6 A (n,n A ) = (u | F) L 2 for all u G H l (R N ). 

Taking u = Ii\ and using the Cauchy-Schwarz inequality, we see that (llOh is satisfied by n A . 
Hence (Vn A ) A>0 is bounded in L 2 and there exist some Q G (L 2 {U. N )) N and a sequence (A n ) ng N 
converging to 0, such that Vn An — 1 Q weakly in L 2 . Note that this implies that Q satisfies 
div (aQ) = div F in the distributional sense, and also that Q is the gradient of some tempered 
distribution n. Besides, we have 

||Vn|| L 2 = \\Q\\ L 2 < liminf ||Vn An || L2 < aZ l \\F\\ L 2. 

As regards uniqueness, it suffices to check that the constant functions are the only tempered 
solutions with gradient in L 2 which satisfy ([8]) with F = 0. So let us consider n G S' with 
Vn G L 2 and div (oVII) = 0. We thus have 

(11) JaVu-VUdx = for all u G H 1 . 

By taking advantage of the Fourier transform and of Parseval equality, it is easy to check that 
for n > 0, the tempered distribution U n := (Id — x( n D))H (where the cut-off function \ nas 
been defined in Section [2]) belongs to H . Hence one may take u = H n in (lllh and we get 

J aVn • Vn n dx = for all n > 0. 

As Vn n tends to Vn in L? and a > a* > 0, this readily implies that Vn = 0. □ 

Let us now establish higher order estimates. 

Proposition 7. Let 1 < p < oo and 1 < r < oo. Let a be a bounded function satisfying ([9]) 
and such that Da G Bp~ l for some s > 1 + N/p or s > 1 + N/p if r = 1. 



s-l 

p,r 

If 1 < p < oo, a G (1, s] and VU G B^ r satisfies ([8]) for some function F such that 
div F G Bp^ 1 then we have for some constant C depending only on s,a,p,N, 



o*||Vn||fl S _ r < C( ||divF|| B a-i +a*(l + a^WDaWzs-i) ||Vn Ni p 



If 2 < p < oo and F is in L 2 and satisfies div F G B^ 1 for some a G (1+N/p— N/2, s] 



then Equation ([8|) has a unique solution H (up to constant functions) such that Vn G 
L 2 D Bp r . Furthermore, Inequality (|10p zs satisfied and there exists a positive exponent 
7 depending only on a, p, N and a positive constant C depending only on s,a,p,N 
such that 

a*||vn|| BpV < cfwdivFWje-i + (l + o^llDoll^) 7 !!^). 
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• If a > 1 and 1 < p < oo then the following inequality holds: 

aJ|Vn||£ CT < C( lldivFll RCT -i + ||Va||z/x> II Villi R <r-i + ||VII||i<»||Vo|| R o--i 

p,r V 1 ' fp,r -Op,r -£3p,r 

Proof. Throughout, (c g ) g >_i denotes a sequence in the unit sphere of £ T . 
The proof relies on two ingredients: 

(i) the following commutator estimates (see Lemmas [6] and [7] in the appendix) 

(12) ||div[o, A 9 ]Vn|| LP < Cc^-^-^HVall^-iUVnll^-i, 

(13) ||div[a, A g ]Vn|| LP < Ccg2 _9(<7_1) (|| Vn|| L oo || Va|| B a-a + ||Vo||£«.||Vn|| B a-i) 

which hold true whenever a £ (0, s] and (s,p, r) satisfies Condition (C) (as regards 
([12]) ) and whenever <r > 1 (as concerns (fl3j) ): 
(iij a Bernstein type inequality (see Lemma [8] in the appendix). 

For proving the first part of the lemma, apply the spectral cut-off operator A 9 to (jBJ). We get 

-div (aA 9 Vn) = div A q F + div ([A„ a] VII) for all q > 0. 
Hence, multiplying both sides by |A g n| p_2 A 9 n and integrating over WL N , we get 

-J |A,n| p - 2 A g ndiv (aA q VU)dx = J \A g U\ p - 2 A g Udw A q F dx 

+ J |A (? n|^ 2 A (? ndiv([A 9 ,a]Vn) dx. 

Apply Lemma [5] to bound by below the left-hand side of the above inequality. Using Holder's 
inequality to handle the right-hand side, we get for all q > 0, 

(14) aa 2q \\A q U\\ p LP < qiAJIII^YlldivAgFllzP + ||div [A q , a]VU\\ LP ) . 



To deal with the last term, one may now take advantage of Inequality (]12p . Since, for q > 0, we 
have ||AqVn||iP 2 9 || A g n||ip according to Lemma [U we get after our multiplying Inequality 
CH) by 2 q ^-^ : 

a*2 9CT ||A 9 Vn|| LP < C7(2 9(CT - 1) ||A (; divF|| L p + c g ||Va|| Brr i||Vn|| B a 7 i) for all q G N. 

Taking the t r norm of both sides and adding up the low frequency block pertaining to A_iVn, 
we get 

(15) aJ|Vn|| B - < C( ||divF|| R a-i + II Vail R ,-i ||Vn|| RCT -i + aJ\ A_iVn|| LP V 

Observe that ||A_iVn||£j> < C||Vn||/,p, and that the following interpolation inequality is avail- 
able (recall that < a — 1 ) : 

livnii^-x <c||vn||i||vn|| B .;. 

Then, applying a suitable Young inequality completes the proof of the first part of the proposi- 
tion. 

Let us now tackle the proof of the second part of the Proposition. As F E L 2 , the existence 

3 a 

P,r 



of a solution Vn in I? is ensured by Lemma [2j Let us admit for a while that Vn G B° and 



let us prove the desired inequality. As p > 2, we have 

T Z j. X P 2' 

^ ' J^p, CO 

Hence, as B°~ 1 is an interpolation space between B Pt0 £ 2 and B° r (here comes the assumption 
that a — 1 > N/p — N/2), one may write for some convenient exponent 6 = 9(p, a, N) £ (0, 1), 

HK7TTM , < n\\\7Tf\\ e . Jivnii 1 - 6 
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In addition, as p > 2, Bernstein's inequality implies that 

||A_iVn|| LP < c||vn|| L 2. 

Hence, plugging the last two inequalities in (fT5|) and using (fTUj) yields 

a*||Vn|| B . r < c([|divF||^- 1 + \\F\\ L 2 +a- 1 ||F||« 2 ||Va|| Brr i(a,||Vn|| B . r ) 1 - e ). 

Then applying Young's inequality completes the proof. 

Remark that Inequality (j!5)l remains valid whenever VII is in B°~ \ Starting from the fact 

iV(i-i) 

that the constructed solution VII is in B p ,o£ , a straightforward induction argument allows 
to state that VII is indeed in B° r . This completes the second part of the proof. 

For proving the last part of the proposition, the starting point is Inequality (|14p which implies 
that 

a*2 9,T ||VA 9 n|| LP < C^^dlAgdivFlliP + ||div [A„ a]VII|| LP ). 

Now, taking advantage of Inequality (|13|) then summing up over q > —1, we readily obtain the 
desired result. □ 



4. Proof of the first local well-posedness result 

As a preliminary step, let us observe that System (fTJ) is time reversible. That is, changing 
(t, x) in (— t, —x) restricts the study of the Cauchy problem to the evolution for positive times. 
To simplify the presentation, we shall thus concentrate from now on to the unique solvability of 
the system for positive times only. 

In the first part of this section, we establish the uniqueness part of Theorem [TJ When proving 
existence, it is convenient to treat the two cases p > 2 and p < 2 separately. The reason why 
is that the proof strongly relies on Proposition [7] which enables to compute the pressure only if 
p>2. Indeed, if p < 2 then only an a priori estimate is stated. 

So, in the second part of this section, we prove the existence in the case p > 2. The third 
subsection is devoted to the proof of Theorem [2] in the case p > 2. It will be needed for proving 
the existence part of Theorem [1] in the case p < 2. The following part of this section is devoted 
to the proof of Theorems [T] and [2] in the case p < 2. In the last paragraph, we justify the claim 
pertaining to the case p > N/(N — 1) (see just after the statement of Theorem [1]). 

For expository purpose, we shall assume in this section and in the rest of the paper that 
r < oo. For treating the case r = oo, it is only a matter of replacing the strong topology by 
weak topology whenever regularity up to index s is involved. 



4.1. Uniqueness. Uniqueness in Theorems [U is a consequence of the following general stability 
result for solutions to ([I]). 

Proposition 8. Let VIIi) and (^2)^2, VII2) satisfy ^ with exterior forces f\ and fi- 

Assume in addition that p\ and P2 are bounded and bounded away from zero, that 5a := u<i — u\ 
and 5p := P2 — p\ belong to C 1 ([0, T\\ L 2 ), that 5f := /2 — f\ is in C([0, T];L 2 ) and that VIIi, 
Vpi and Vui belong to L 1 ([0, T]; Then for all t £ [0, T], we have 

L2 + ||(^4i)(*)ll^ < e A{t) (\\mh* + IKvWCO)!^ + [ e- A ^\\^p 2 Sf\\ L2 



(16) 

with A(t) : 







V^2 



+ 



P1VP2 



+ || Vni ll^oo dr 



Proof. On the one hand, as 
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taking the L 2 inner product with 5p and integrating by parts in the second term of the left-hand 
side yields 



V/91 



(17) WWlli? < \M0)\\u + / \\(Vp 2 Su)\\l* - 

Jo \fP 2 

On the other hand, denoting V5Q := VII2 — VLTi, we notice that 



dr. 

L°° 



p 2 (d t 5a + u 2 ■ V6u) + V&I = p 2 [5f + — — Vn x - 5u ■ Viti 

V P1P2 

So taking the L 2 inner product of the second equation with Su, integrating by parts and using 
the fact that div 5a = and that 

d t p2 + u 2 ■ V/9 2 = 0, 

we eventually get 

\\(yfpM(t)h* < \\(VPoSuMh*+ f%\Vp 2 Sfh*MMv T +\\Vu l \\ L ^\\^p 2 5a\\Adr. 



L° 



Adding up Inequality (|17p to the above inequality and applying Gronwall lemma completes the 
proof of the proposition. □ 

Proof of uniqueness in TheoremUl Consider two solutions (pi,ui,Vrii) and (p2;W2iVIl2) 
of dH) with the same data. Under the assumptions of Theorem [lj it is clear that the ve- 
locity and pressure fields satisfy the assumptions of the above proposition. As concerns the 
density, we notice that Ui G C([0,T]; L 2 ) and V pi G C([0, T];L°°) for i = 1,2 implies that 
d t pi £ C([0,T];L 2 ). Hence we have 5p G C 1 ([0, T];L 2 ). Therefore Inequality CE]) implies that 
(pi,ui,VLIi) = (p 2 ,u 2 ,VU 2 ) on [0, T] x R-^. This completes the proof of the uniqueness in 
Theorem [TJ 

4.2. The proof of existence in Theorem [1} the case 2 < p < 00. We notice that, formally, 
the density-dependent incompressible Euler equations are equivalent tcEl 

{dta + u ■ Va = with a := 1/p, 
d t u + u-Vu + aVn = /, 
-div (aVn) = div (u ■ VVu) - div /. 
Let us give conditions under which this equivalence is rigorous. 

Lemma 3. Let u be a time- dependent vector-field with coefficients in C 1 ([0,T] x WL N ) and 
such that Qu G C 1 ([0, T]; L 2 ). Assume that VII G C([0, T];L 2 ). Let p be a continuous bounded 
function on [0, T] x K N which is positive and bounded away from 0. 

// in addition divn(0, •) = in H N then (p,u,VU) is a solution to ([1]) if and only if 
(a, u, VIT) satisfies (fl~8j) . 



Proof. If (/?, u, VIT) satisfies ([T]) then, owing to p > 0, we see that a := 1/p satisfies the first 
equation of (fTSj) . Next, applying Operator div to the velocity equation of ([T]) divided by p, 
and using that Vu = u yields the third equation of (fTHj) . 

Conversely, if (a,u, VII) satisfies (|18p . it is obvious, owing to positivity, that p := 1/a satisfies 
the density equation of ([I]). In order to justify that the other two equations are satisfied, it is 
only a matter of proving that divn = 0. For that, one may apply Q to the second equation. 
Then, using the third equation, we discover that 

d t Qu + Q(u-VQu) = 0. 
Recall that Qu G C 1 ([0, T];L 2 ). Therefore, taking the L? inner product with Qu, we get 

\j t \\Q u \\h + (Q( u ■ v 2^) I Q") L 2 = °- 
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As Q T = Q and Q 2 = Q, we thus get after integrating by parts in the second term: 

d ..o f ,2 

^■||Qm|| L 2 = / \Qu\ div udx, 
and, as Qu(0, •) = 0, Gronwall lemma ensures that Qu = 0. Hence divu = 0. □ 

As explained in Lemma El it suffices to solve System (fT8|h So, for T > 0, let us introduce the 
set Et of functions (a, u, VII) such that 

a G C 6 ([0,T] x R N ), Va G C([0, T]; B'- 1 ), 

u G ^([0, T]; L 2 ) n C([0, T]; B£ r ), VII G C([0, T]; L 2 ) n tfflO, T]; B£ r ). 

We denote 

a* := inf ao(x), a* := sup ao(x), := inf /9o(s) and p* := sup /?o(^)- 

xGM^ zeRlv xGM^ xeR iV 

Note that if p is bounded and bounded away from zero, and satisfies Vp G then the same 

properties hold for a (and conversely). This may be easily shown by combining Propositions 
[2] and [H Moreover, there exists some constant C depending only on a*, a*, N and on the 
regularity parameters such that 

C~^Vp\\ B s-i<\\Va\\ B s-i<C\\V P \\ B s ? i. 
This fact will be used repeatedly in the rest of the paper. 

Step 1. Construction of a sequence of approximate solutions. As a first step for solving (fT8|) . we 
construct a sequence (a n , u n , VII n ) ng N of global approximate solutions which belong to Et for 
all T > 0. 

For doing so, one may argue by induction. We first set (a , u°, VIT ) := (ao,^o,0). Next, 
we assume that (a n ,u n ,VII n ) has been constructed over M + , belongs to the space Et for all 
T > and that there exists a positive time T* such that for all t G [0,T*], 



(19) a* < a n (t,x) < a*, 

(20) ||Va n mi| RS -i < 2||Va ]| RS -i for all t G [0, T*l, 

(21) ||y^i«(i)|| i2 < v^(4||V^uo|U 2 +8^11/11^^)) with p":=l/a", 

(22) U n (t) < 4U (t) + Cop^olldiv/H^-i) + C (p*A r +1 (\\u \\ L 2 + ||/|| L i (L2) ), 

«*l|Vn«|| Ltl( ^ r) <c(| o \^(r)) 2 dr + ||div/|| itl(Brrl) 

(23) +(P*A y(\\u \\ L 2 + \\f\\ Ll{L2) )y 

(24) [|VIP||ii (La) oIt < ^||y^« || L 2 + 



with A := a* + \\Da \\ B s-x, U (t) := \\u \\ B ^ + II/IIl^,,) an d t/ n (t) := The 
positive exponent 7 is given by Proposition [7J The constants Co and C depend only on 
(s,p, r) and N, and may be made explicit from the following computations (in fact one can take 
Co = 2C 2 with C large enough). 

Denoting by ^ n the flow of u n , (which belongs to C 1 (M+ x R N ) owing to u n G C(R + ;B^ r ) 
and to Bp~ l ^ C^), we set 

a n+1 (t,x):=a ((W)- 1 (x)) and p" +1 (i, s) := po(W4TV))- 
As V'f is a diffeomorphism over for all t > 0, we have 

11 /•j.mi 11 _ 11 „* j II _n+l/.iMI ll„ II „* 
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Hence (fT9j) is satisfied by o n+1 . In addition, we have 

d t a n+1 + u n • Va n+1 = 

so that for all i € {1, • • • , TV}, 

<9id;a n+1 + u n ■ Vd ia n+1 = -d iU n ■ Va n+1 . 

As dia? t ~t = dido G -Bp^ 1 by assumption, (a slight generalization of) Proposition [5] combined 
with Gronwall lemma guarantees that Va n+1 E C(M + ; Bi^) and that 

(25) HVa^WH^ < e^M^IV^x. 
So if we assume that T* has been chosen so that 

(26) C [ T U n {t)dt < log 2 

J o 

then a n+1 satisfies (j20l) . 

Next, we want to define u n+1 as the unique solution in C(R + ; B s p r ) of the transport equation: 

(27) 8 t u n+l + u n ■ Vu n+1 = -a n+1 Vn n + /, upj = u . 

That the right-hand side belongs to Lj oc (M + ; Bp r ) is a consequence of Corollary Q] and of the 
embedding B^} > L°°. In addition, we have for a.e. positive time 

(28) |K +1 VlT||^ r < C(\\a n+l \\ Laa + ||Va ra+1 || B ^i)||Vn™||^ r . 

So finally, the existence of u n+l £ C(R + ; r ) is ensured by Proposition [5] and we have 

II n+lf,\\\ ^ C (lu n (r)d,T ( n I, 



(29) + /^^^^(^^'((lla^lUco + ||Va" +1 || B| -i)||Vnl B|r + dr 

Therefore, if we restrict our attention to those t that are in [0, T*] with T* satisfying (|26p . we 
see that for all t€ [0,T*] 5 

< 2t/b(*) + CA / ||Vn n || B| r dr with A := a* + ||Va || B «-i. 

« ' 

So if we assume that T* and Cq have been chosen so that 

pT* 

(30) 2C 2 p*A U n (t)dt<l andC = 2C 2 

JO 

then taking advantage of Inequalities ([23]) and ([22]) , we see that ii n+1 satisfies ([22]) on [0,T*]. 

Let us now prove (|2ip for -u n+1 . First, we notice that the right-hand side of (|27p belongs to 

C(M + ;L 2 ) so that u n+l is in (^(R+jZ 2 ). As p n+1 is bounded and C 1 with respect to the time 

and space variables, this allows us to take the L 2 inner product of the equation for u n+l with 
pn+lyu+i _ ^ reac jiiy g e t 

(31) -^-\\^/p~^ u n+1 \\ 2 L2 - f p n+1 u n+l -fdx = - I p n+1 \u n+1 \ 2 divu n dx - (Vn n | u n+l ) L2 . 
2 ctt J 2 J 

Let us point out that u n and u n+l need not be divergence-free, so that the right-hand side may 
be nonzero. However, from the above inequality, it is easy to get 

\\{47^u n+i ){t)\\ L ,<\\^- P ,u Q \\ L , 

+ /Yv^iViT 1 !^ + x/pII/II^ + i|K/^^ +1 |L2||div^|| L oc') dx. 
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So, if we assume that C has been taken large enough in ([26]) then Gronwall's lemma implies 
that 

(32) \\{J^u n+l ){t)\\ L * < 2(|| Vpo^oII^ + v^ll/ll^ 2 ) + v^iViF||^ 2) ). 

Now, putting the above inequality together with Inequality (|24p ensures that Inequality (|2ip is 
also satisfied by u n+1 on [0, T*]. 

To finish with, we have to construct the approximate pressure II n+1 . For that, we aim at 
solving the following elliptic equation 

(33) div (a n+1 VII n+1 ) = div (/ - u n+l ■ VVu n+1 ) 
for every positive time. 

We have already proved that a n+1 satisfies the required ellipticity condition through (|19p . 
Moreover, as u n+1 G C(R + ;B°_ r ), Remark [2] ensures that VVu n+1 is in C(R + ;B^~ 1 ). As 
B s p - 1 ^ L°° and u n+1 G C(R + ; L 2 ), we thus have u n+l ■ VVu n+1 G C(R+; L 2 ) and 

\\u n+1 • VVu n+1 \\ L 2 < v ^||Vi^+ T « n+1 || L 2||VW +1 || J Lo , 
< CV^|| v / p s + T M n+1 || x2 ||Vu n+1 || B «-i. 

Therefore Lemma [2] guarantees that (|33|) has a solution VII n+1 in C(R + ;L 2 ) which satisfies 

(34) ^||Vn" +1 || Ltl{L2) < ||/|Li (L2) + Cv^ fu n+1 \\^^u n ^\\ L ,dr. 

J 

Let us insert Inequality (f3"2"j) in the above inequality. We see that if T* has been chosen so that 

(35) ACa*p* [ T U n+1 dr < 1 

Jo 

then Inequality ([34"]) implies that 

l|vn™ +1 || itl(i2) < l P *\\f\\ Ll(L2) + IVfWVtovoWv + ^livn n || itl(i2) , 

hence Inequality dMD is satisfied by VU n+1 on [0,T*]. 

In order to prove that VIl n+1 belongs to Lj oc (R + ;Bp r ), one may apply the second part of 
Proposition [7J Indeed, because, owing to div Vu n+1 = 0, we have 

div (u n+1 ■ VVu n+1 ) = Vu n+1 : VVu n+1 
and as B^} is an algebra, the term div (u n+1 • Wu n+1 ) is in Btlj, and 

||div(u n+1 • VVu n+1 )\\ B s- } < C{U n+1 f. 
Hence Proposition [7J implies that for all t G R + , 

a*l|vn^|| Ll{B|ir) < c(J\u^) 2 dT + ||div/|| itl(B|rrl) 

+ (1 + p* llDa^X^^ym^ + ||u^ • VPn" +1 || itl(L2) )), 
whence, using ([20]) at rank n + 1 and Holder inequality, we get 

a*||Vn^|| itl( ^ r) < c(j\u^ l ) 2 dr + [Idiv/H^-ij 

+ (p^o) 7 (||/||L t ML 2 ) + v^||vWn" +1 || Lr(L2) jTt^dr)). 
Taking advantage of Inequality ([2ip at rank n + 1 one can now conclude that if (|35[) holds then 
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At this stage we have proved that if Inequalities (|19l) to (I24p hold for (a", u n , VIP) then they 
also hold for (a n+1 ,u n+1 ,VU n+1 ) provided T* satisfies Inequalities fl2(|, ([30]) and (J35]). Note 
that pop is the strongest condition. Obviously it is satisfied if we set 
(36) 

T* := sup{t > 0/p*tA o (lI (t) + p*A \\dwf\\ Ll{B s- 1) + (p*A y +1 (\\u \\ L 2 + ||/|| L i (L2) )) < c} 
for a small enough constant c depending only on s, p and N. 

Step 2 . Convergence of the sequence. Let a n := a n — ao- In this step, we shall establish that 
(o n ,/,Vn n ) na is a Cauchy sequence in C([0, T*];L 2 ). 

Let 6a n := a n+1 - a n , 5u n := u n+1 - u n and CP := n n+1 - IF. We have for n > 2, 

!a i( ii n + u 11 ■ V6a n = -5a n ~ l ■ Va n , 
d t 5u n + u n ■ V6u n = -A/ 1 " 1 • Vu n - a n VdI n - 1 - 6a n VU n , 
div (a n ' 1 V6Tl n - 1 ) = -div (8a n ~ l ■ Wu n + u n ~ l ■ VVSa 11 ' 1 + fo^VIF). 

For all n E N, we have dtd n+l = —u n ■ Va n+1 . So, given that, according to the previous step, 
u n E C([0,T*];L 2 ) and Va n+1 E C 6 ([0,T*] x R N ), and that aj£J = 0, we discover that a n+1 , 

and thus also 6a n , are in C 1 ([0, T*]; L 2 ). Taking the L 2 inner product of the equation for 5a n 
with 6a n , we thus get 

' ,! I / '(&™) 2 div u n dx- J du^ 1 -Va n &L n dx, 



\& n \\h 



2dt" " h 2 
whence for all t £ [0,T*], 

(38) \\Sa n (t)\\ L 2 < ~ r||divn n || L oo||& n || i2 (ir+ t \\V a n \\ Laa \\5u n ~ l \\ L 2 dr. 

2 Jo Jo 

Next, taking the I? inner product of the equation for 5a n with p n+1 5u n , performing integration 
by parts and using the equation for p n+1 , we get 

i|/,*>»l^ = i/ p « + > n ^v»»^-/,«&».(&»-'w + a«v^-' + &»vn»)^. 

Hence 

\\^p^\t)du n (t)\\ L2 < \ [\\&ivu n \\ L ™\\^fp~^5a n \\ L ,dT 



2 

t 

n— in - i II -nil. llrrJB-m-ll 



(39) + / (||Vu n || LO o|| ( Si n - 1 || L2 + ||a n || L <»||V«5n n - i || L 2 + ||VlF|| L oo ||& n || L2 ) dr. 

Jo 

Adding up Inequalities (|38p and (|39p . applying Gronwall lemma and using the fact that p n+1 > 
p* and the bounds stated in the first step, we thus get for all t E [0, T*], 

(40) \\(5a n ,6a n )(t)\\ L 2 < C T J f \\{5a n -\5a n - x ){r)\\ L , dr + f \\VSi n -\r)\\ L 2 dr 

\Jo Jo J 

where the constant Ct* depends only on T* and on the initial data. 

In order to bound V<Si n_1 , we shall use that for any C 1 vector-fields a and b, we have 
div (a • V6) = div (b ■ Va) + div (a div b) — div (6 div a). 

Applying this to a = u n ~ 1 and b = V5u n ~ l and bearing in mind that &\~vV5a n ~ l = 0, we 
deduce from the third equation of (|37p that 

div (a n - l VM n - 1 ) = div (7'<5u*- 1 divti n - 1 - V5a n - 1 ■ Vu n ~ l - da 71 ' 1 ■ VVu n - Sa n ^VU n ). 

Therefore, Lemma [2] and the fact that HT'H/^l 2 ;!, 2 ) = 1 guarantee that 

a^llV&I™- 1 !!^ < [|& i n - 1 || L a(||div« ,l - 1 ||r<» + llVn"- 1 ^^ + ||VPu n || L oo) + H&™- 1 || £2 ||Vn n || Loo . 
Using the uniform bounds of the previous step, we thus get for all t E [0, T*], 

r ai\ llvTjTTn-ln ^ r< l\\ SL.n— 1 II i ll.cn— 111 \ 
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Plugging Inequality (HTj) in Inequality (HOI) , we end up with (up to a change of Ct* ) , 

\\(&L n ,du n )(t)\\ L 2 <C T * f'w^-^du^^W^dr. 

Jo 

Arguing by induction, one may conclude that 

sup \\(Sa n ,6u n )(t)\\L* < { T , ] sup ||(& >°)(t)|| L2 . 
te[o,T*] n! te[o,T*] 

It is now obvious that both (a n ) nG N and (tt" - ),,,^ are Cauchy sequences in C([0, T*]; L 2 ), hence 
converge to some functions a and u in C([0,T*]; L 2 ). Taking advantage of (f4Tj) . it is also clear 
that (Vn n )„gN converges to some function VII in C([0, T*}; L 2 ). 

Step 3. Final checking. Let a := ao + a. We now have to check that (a,u, VII) is indeed a 
solution to (HJ) and that it has the properties stated in Theorem [TJ From the previous step, we 
already know that (a — ao), u and VII are in C([0,T*];L 2 ). Moreover: 

• As (Va")„ g fj is bounded in L°°([0, T*]; B*~ ) and as Besov spaces have the Fatou prop- 
erty, we deduce that Va belongs to L°°([0, T*]; Bp" 1 ). Since (a n ) ng N is bounded in 
L°°([0,T*} x R N ), we also have a £ L°°{[0,T*} x R^). 

• As (n n ) neN is bounded in L°°([0, T*};B^ r ), we deduce that u £ L°°([0, T*];B^ r ). 

• Finally, as (VII n ) rag N is bounded in -L 1 ([0, T*]; B!l r ) we deduce that VII belongs to 
H([0,T*];B' r ). 

Arguing by interpolation, we see that the above sequences converge strongly in every interme- 
diate space between C([0, T*]; 1?) and C([0, T*]; ) which is more than enough to pass to the 
limit in the equations satisfied by (a n ,u n ,VII n ). Hence (a, u, VII) satisfies (I18p , 

Passing to the limit in ([3Tj) . we see that, in addition, (p,u) satisfies the energy equality (J3j) - 
Finally, the continuity properties of the solution with respect to the time may be recovered 
by using the equations satisfied by (a, u, VII), and Proposition [5j 

4.3. A continuation criterion. The key to the proof of Theorem [2] is the following lemma: 

Lemma 4. Let (s,p,r) satisfy Condition (C) with 1 <p < oo. Consider a solution (p, u, VII) 
to ([I]) on [0,T[xR N such thaE u G C([0,T); B^ r ) and 

P*<P<P*, peC([0,T) xR N ) and V p £ C([0,T); B^ 1 ). 

If in addition 




Proof. Note that a := 1/p satisfies the same assumptions as p. Therefore we shall rather work 
with a, for convenience. Recall that 

(43) dtdka + u ■ V<9fca = — d^u ■ Va for k = 1, • • • , N. 

So, applying Operator A q to the above equality and using that div-u = 0, one may write (with 
the summation convention) 

dfAqdka + u ■ VA q dka = —A q (dkU ■ Va) + dj[u J , Agjc^a. 

Therefore for all t € [0, T), 

(44) \\A q d k a(t)\\ LP < \\A q d k ao\\ LP + / || A q (d k u ■ Va)\\ LP dr + / \\dj[u j , A q ]d k a\\ LP dr. 

Jo Jo 
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According to Proposition El the term d k u • Va belongs to Bp J. 1 and satisfies 

\\dkU ■ Va\\ B s-i < C(\\dku\\Loo ||Va|| BS -i + ||Va||/,oo ||9fcii|| BS -i) 
while Lemma [7] ensures that for all q > —1, 

\\dj[u j , A q }d k a\\ LP < Cc q 2 q ^ 1 \(\\d k a\\ L ^\\Vu\\ B s-i + ||Vu||£°o ||<9 fc a|| B s-i). 

Using the definition of the norm in -B*" 1 , we thus get after summation in (|44|) that 

(45) \\Va(t)\\ B s-x < [IVooIIb.-! +Cjr t (||V«|Ua ||Vo|| B .-i + ||Va|| z =o||V«||^-i) dr. 

In order to bound the velocity, let us apply the last part of Proposition [5] to the velocity equation, 
and the following inequality (which stems from Corollary [T|) : 

||aVn|| B s r < C7(a*||Vn|| B s r + ||VII||i<x.||Va|| BS -i). 
We get for all t G [0, T), 

\\umB^ r <e C Io^- d ^\ Uo \\ B , r 

(46) + jf* e~ c Io ll v «H^ ^' (||/|| fl . r + Ca*\\ Vn|| B . r + C||Vn|| L . ||Va|| B| -x) dr) . 

In order to bound the pressure term, one may use the fact that 

div (aVLT) = div / — div (u ■ Vu) 
and apply the last part of Proposition [7J Performing a time integration and using the fact that 
|| div (u ■ X7u)\\ B s-i = || Vu : Vtt[| B »-i < C|| V«[|i,«> || Vu|| B s-i , 

we get 

a*l|Vn|| L i (B|r) < C^lldiv/H^^-i) 

+ / (||Vii||L°° II Vtill R »-i + ||Va||ioo ||VII|| R «-i + || Vn||i°o || Vo|| R «-i) dr. 

Let us insert this latter inequality in (|46p • Then adding up Inequality (|45p and applying Gronwall 
lemma we end up with 

||Va(i)|| B |-i + ||tt(t)||B», r < Cex Af Q ||(Va, V«, Vn)||ioc dr^j 

(47) (l|Va || B |-i + ||«o||b- p + \\fhUBi, r ) + f* l|Va|| L oc||Vn|| Brrl dr) 

for some constant C depending only on the regularity parameters and on N, a* and a*. 

Now, let us notice that Va is bounded on [0, T) x R . Indeed, from Equation ([13]) and 
Gronwall lemma, we see that 

||Va(t)IU~ < e /oll v «ll^||Va ||L~. 

As Vn is in Z^flO, T); B^ 1 ) and Vu is in L 1 ([0, T); L°°) by assumption and as B 3 ^ 1 ^ L°°, 
we discover that both the last term in (|47p and the exponential term are bounded on [0, T). 
This completes the proof of the lemma. □ 

The following lemma implies the first part of Theorem [2] in the case p>2. 

Lemma 5. Let (s,p,r) satisfy Condition (C) with 2 <p < do. Consider a solution (p, u, VII) 
to © on [O^fxM^ suc/i i/iafl 

• P*<P<P*, p £ C([0,T) x R N ) and Vp£ C([0,T); B*' 1 ), 
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• neC([0,T);^ r )nC 1 ([0,T];L 2 ), 
. VnGC([0,T);L 2 )nL 1 ([0,T); J B^). 

If in addition Condition (I42p is satisfied then (p, u, VII) may be continued beyond T into a 
solution of ([!]) with the above regularity. 

Proof. Lemma H] ensures that |M|i,°°(.Bs p ) and \\Da\\ are finite. So one may set 

e := c(p*A r 1 (u (T) + p*A \\div (B .- pl) + (p*A y +1 (||« || i2 + ||/||^ (I/2) )) _1 

where c is the small constant (depending only on N and (s,p,r)) defined in (|36p . 

Then we know from the proof of Theorem [1] in the case p > 2 that for any T 1 < T, System (pQ) 
with data (p(T'), ti(T'), f(T' + •)) has a unique solution up to time e. Taking T' = T — e/2 we 
thus get a continuation of (p, u, VII) up to time T + e/2. □ 

Let us now justify the last part of Theorem [21 It stems from the following logarithmic 
interpolation inequality (see e.g. [15J): 

IIVuIIl 00 < C{ 1 + II Vnll ro log(e + II Vu|| ns-i) ) with ||Vu||oo := sup ||w(2 _i? L')Vm||loo 

V -°oo,oo v -°p,r 'J -°oo,oo QgjZ 

which holds true whenever the embedding of B^ 1 is not critical (that is s > 1 + N/p). 

Then, arguing exactly as in Proposition 5.3 of [TT], we discover that Condition (|42p may be 
replaced by the following weaker condition: 

r T 

(48) / (llVitlloo + II Vnll o—i) dt < oo. 

Now, it is classical (see e.g. [2], Chap. 7) that there exists some constant C such that 

||Vu||do < C||curln||Loo. 
This completes the proof of Theorem [2] in the case p > 2. 

4.4. The case 1 < p < 2. Note that by virtue of Proposition [2l the data satisfy the assumptions 
of the theorem for the triplet (s — N/p + N/2,2,r). Hence, applying the theorem in the case 
p = 2 supplies a local solution with the B2 r N / p+N l 2 regularity. However, proving that the Bt, r 
regularity is also preserved, is not utterly obvious. For proving that, we shall proceed as follows: 

i) first, we smooth out the data so as to get a solution in H°° := (~) a H a for which the Bp r 
regularity is also preserved; 

ii) second, we establish uniform bounds in Bp r on a fixed suitably small time interval; 
Hi) third, we show the convergence of the sequence of smooth solutions and that the limit 

has the required properties. 

Step 1: smooth solutions. Set 

Oq := S n a , Uq := S n u and f n := S n f 

where S n is the low frequency cut-off introduced in Section [2 
Note that for all large enough n G N and t £ [0, To], we have 



(49) a*/2 < a% < 2a*, \\Da%\\ B s-i < C\\Da \\ B s-i, 

(50) II«oIIb» p < c\\u \\ B s tr , 

(51) \\nLUB°, r )<C\\fhUB kr ) and \\nt)\\ L *<\\f(t)\\ L 2. 

It is also clear that (with obvious notation) Voq and Uq are in B2? r (hence also in H°°) and 

a.i m s- nl\c\ rn 1. uoo\ r\ t 1/rn ml. doo \ 
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Finally, taking advantage of Lebesgue's dominated convergence theorem one may prove, if 
r < oo, that@ 

Do,q — > Dao in Bp~ l and (aft — a) — > (ao — a) in L p * with a := 1/p, 
Uq — ► no in B*~ x , 

f n in Li([o,r ];B' r )nC([0,r ];L 2 ). 

As usual, the strong convergence has to be replaced by the weak convergence if r = oo. 

Applying Theorem [1] in the case p = 2 and using the fact that the lifespan does not depend 
on the index of regularity (see Remark [I]), we get a local maximal solution (a n ,n n , Vn n ) with 
Da", u n and VIP in C([Q,T*);H°°), and 

(52) a*/2 <a n < 2a*. 

Note that as a n and p n are just transported by the (smooth) flow of u n , we also have 

(53) \\(p n (t)-p)\\ LP * =\\Po -p\\lp' <C\\ P o-p\\lv* for all ie[0,T*) 

(and similarly for a n ) and Va n and V p n belong to C([0,T*); B^ 1 ). 

Let us now establish that VIP 1 is in L 1 ([0, T];B^ r ) for all T G [0,T*). Fix some T G [0,T*). 
Applying Operator div to the momentum equation of (pQ) and using that div u n = yields 

(54) AIT = div (p n f n ) - div (p n u n • Vu n ) - Vp n ■ d t u n . 

According to Proposition [H F i— » D 2 U is a self-map on Bp r . Hence, in order to show that 
Vn n G ^([CT];^), it suffices to establish that Vn n e L l ([0, T]; D>) and that An™ G 
L\%T]-B s - r x ). 

Let us first show that all the terms of the right-hand side of (|54p are in L ([0, T]; Bp~ ). 
Since, by assumption, / n G L ([0, T];!?*) and as it as been established that p n G L°° and 
Vp n G -Bp"; 1 , Corollary Q] implies that div (p n f n ) G L x ([0,T]; -B^ 1 ). For the next term, we use 
that for all i G {1, • • • , N}, 

{p n u n ■ Vu n T = J2 T' pn{un)j d jU n + T d]{uny p n {u n y. 
j 

By embedding, p n u n and Vu n are in LP (recall that p* > 2 > p) and, arguing as for p n f n , 
one can check that p n u n is in H°°. Of course, Vn n is also in H°°. Given that 1/p = 1/p* + 1/2, 
continuity results for the paraproduct and remainder in the spirit of Proposition [3] (see |18j ) 
ensure that p n u n ■ Vu n is in Bp~ . 

For the last term in (|54p . one may write that 

Vp n ■ d t u n = T' dtun ■ Vp n + T Vp n • d t u n . 

As, by embedding, dtu n G and as Vp n G Bp~ x , continuity results for the paraproduct 
ensure that the first term in the right-hand side is in Bp/J, 1 . Concerning the second term, one 
may use that Vp n G LP* (by embedding) and that dtu 11 G H°° (from the equation). Hence 
An n is indeed in L ([0, T];Bp~ ), as claimed above. 

In order to establish that Vn n G L 1 ([0, T];LP), we use the fact that, owing to div dtvP = 0, 
one may write 

(55) An n = div (p n / n - p n u n ■ Vu n - {p n - p)d t u n ) . 

Hence, it suffices to check that p n f n , p n u n -Vu n and (p n —~p)dtu n are in L 1 ([0, T]; L p ). For p n u n 
this is obvious as, by embedding, f n G L 1 ([0, T]; L p ) and u n G C b ([0,T] x R N ). By embedding, 
we also have Vu n G C([0,T];L°°) and u n G C([0,T]; LP), hence ,o n n n • Vu n is in L x ([0, T]; LP). 

To deal with the last term in ([55]) the property that (p n — p) G C([0,T];L p ) comes into 
play. Indeed, from the velocity equation, as the solution is in H°° , one easily gathers that dtu n 
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belongs to C([0, T];L 2 ). Hence Holder's inequality (note that 1/2 + 1/p* = 1/p) ensures that 
(p n -p)d t u n eC([0,T];LP). 

To finish this step, one has to prove that u n is in C([0, T*); Bp T ). In fact, from the product 
laws in Besov spaces and the properties of regularity that have been just established for the 
pressure and the density, we get 

d t u n + u n ■ Vu n = f n - a n VU n G LL([0, T*); 2£ r ). 

As Uq G B*, Proposition [S] ensures that u n G C([0,T*); B* ). 



Step 2: Uniform estimates. Let us remark that, by Sobolev embedding and owing to (pE^) . (f5Uj) . 
(foTj) . one may find some index a > d/2 + 1 such that (Dag) n6 N, (ug)neN and (f n ) n eN are 
bounded in H a ~ x , H a and C([0,T ];L 2 ) n L 1 ([0, Tg]; H a ), respectively. Taking advantage of 
Theorem [T] in the case p = 2 and of the lower bound provided by (|36l) we thus deduce that there 
exists some time T > and some M > such that for all n G N, we have T* > T and 

(56) llValjoo^-i) + |K|| i? o (//CT) + ||Vn n || L ^ CT) < M. 

Of course the energy equality ([3]) is satisfied on [0, T] by any solution (a n , u n ,VII"). Recall that 
in addition, according to the previous step of the proof, f|53j) is satisfied and 

Va n G C([0, T]; B^ 1 ), u n G C([0, T];B^ r ) and Vn n G ^([0, T};B s p r ). 

We claim that, up to a change of T, the norm of the solution may be bounded independently of n 
in the space Ex defined in Subsection 14.21 In all that follows, we denote by Cm a "constant" 
depending only on (s,p, r, N, a*, a*) and on M. 
From Proposition [5j we have 

C f* ||V« n |l s _i dr 

(57) ||Vo»(t)||j» < \\Va%\\ B s e Jo 



and, arguing as for proving Inequality (|28j) . 

C f* ||Vu n |l s _i dr ( 

IK(t)||^ r <e Jo " [\\uIWb*,, 

(58) +j Q e Jo (\\p\\ BSpr + + ||Va"|| B| -x) ||VH"|| B| Jdrj . 

In order to bound Vn n , we apply the first part of Proposition [7] to the following equation: 

div (a n Vn n ) = div (f n - u n • Vu n ) . 

Using the fact that -B*" 1 is an algebra and the relation div {u n ■ X7u n ) = Vu n : Vu", we end up 

with 

(59) 

"*l|vni^ r) <c(\\n L]{B$r) + /jK||^,dr + a,(i+ llDa ^ ) Vn"!^^ 

In order to "close the estimate", we now have to bound VH. n in LP . For that, we apply the 
standard LP elliptic estimates stated in Proposition [6] to ([55]) , and Holder inequality so as to get 

\\vn^ (LP) <c( P *(\\p\\ Lj{LP) + J*\\u n \\ LP *\\Vu n \\ L idT^ + \\p n -p\\ L ° ■ " a - n 

Note that, by Sobolev embedding, we have 

IK|| LP . <c\\u n \\ H «. 

So finally, there exists some constant Cm such that 

IIY7TT"II ._ . <r r., , 



(LP*)\\"tU \\Ll(L 2 ) 



,n\\2 7 i „ f~i I 1 I n t \ D P,r ) 
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Plugging this latter inequality in ([59]) . we thus get 

a.||VII»|Uj (fl . r) < c(\\r\\ Lmr) + jf \\VuX^ dr + a,C M (l + 

It is now easy to conclude this step: denoting 

U n (t):=\\u n (t)\\ B s r and ,4 n (i) :- a' + \\Dn" ( 

and assuming that T <T has been chosen so that 



(60) C / \\Vu n \\ RS -i dr < log 2, 

Jo p ' r 

the above inequalities and (09]), ([50]) . ([51]) imply that for all t £ [0, T] we have 

A n (t) < 2A with ^4 := a* + II-DooIIr— i 

and 

E/ n (t) < 2(?7 (t) + Cp*A J*(\\f\\ B s r + (U n f + C M a*Al\ dr. 

So finally, there exists a nondecreasing function F depending only on the norm of the data and 
such that for all t G [0, T], we have 

U n (t) < 2F(t) + Cp*A f\u n {T)) 2 dr. 

Jo 

Therefore, if in addition 

(61) 2CA [ T U n (r)dT < a, 



then we have U n < 4F on [0, T\. 

By arguing exactly as in the case p > 2, it is easy to see that Condition (JBT]) is satisfied if T 
is small enough (an explicit lower bound may be obtained in terms of the data). So finally, we 
have found a positive time T so that (a n , u n , VIT n ) ng N is bounded in the space Et . 

Step 3: Convergence of the sequence. Let 5a n := u n+1 — u n , 5p n := p n+1 — p n and <MI n = 
n n+i _ n n Applying Inequality (p} to the solutions (p n ,u n ,VU n ) and (p n+1 , n n+1 , Vn n+1 ) 
and using the uniform bounds that have been established in the previous step, and (I52p ensures 
that there exists some M > such that for all t € [0, T] and n £ N, we have 

(62) ||* n (t)|U 2 + II^WIIl^ < Af (||* n (0)|| £2 + ||^ n (0)|| L2 + f Q \\5f n \\ L i dr^j. 

Now, from the definition of ag and the mean value theorem, we get for large enough n, 

\\5p n (0)\\ L 2<C2- n \\Da \\ L 2. 

Similarly, we have 

\\5u n (0)\\ L 2 + WSrhlW < C2-"(||^ ||i 2 + \W\\ L i T w). 

So Inequality ([62]) entails that (p n — / o°) ng N is a Cauchy sequence in C([0,T];L 2 ) and that 
(n n ) ng N is a Cauchy sequence in C([0,T]; L 2 ). Then, using for instance ([55]) . we see that 
(VII n ) ng N is also a Cauchy sequence in C([0, T];L 2 ). 

Finally, from the bounds in large norm that have been stated in the previous step, and the Fa- 
tou property for the Besov space, one may conclude that the limit (a, u, VII) to (a n ,u n , VII n ) ng N 
converges to some solution (pQ) and has the desired properties of regularity. As similar arguments 
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Let us now establish Theorem [2] in the case p < 2. Let (p, u, VII) be a solution with the 
properties described in Theorem [U Note that Lemma H] is also true if p < 2. So the only change 
lies in the proof of Lemma [5] which now uses the (new) lower bound for the lifespan that may be 
obtained from the computations of step 2, instead of (|36p . This gives the first part of Theorem [2] 
As in the case p > 2, the last part of the proof of the theorem is a mere consequence of the 
logarithmic interpolation inequality stated in [15] . 

4.5. Removing the assumptions on the low frequency of the data. As pointed out in 
Section [U in dimension N > 3, the supplementary assumption that (po — p) G L p is not needed 
if p> N/(N -1). 

In order to see that, one may repeat the proof of the theorem in the case 1 < p < 2. As 
before, bounding VII n in L T (L P ) is the main difficulty. For that, one may decompose VII™ 
into two terms VII n and VIE; such that 

AIT? = div (p n f n - p n u n ■ Vu n ) and AIL] = Vp n ■ d t u n . 

On the one hand, as before, one may write that 

||vn n || L p < Cp*(\\r\\ LP + \\u n \\ LP * \\Vu n \\ L 2). 

On the other hand, we have 

Vn n = (-A) _1 V(Vp n • d t u n ). 

Recall that in dimension N > 2, the kernel of Operator (— A)~ X V behaves as |x| 1_Ar . Hence, 
according to the Hardy-Littlewood-Sobolev inequality, if 1/p + l /N < 1 then we have 

||Vn^|| LP < C\\Vp n ■ d t u n \\ Lq with 1 = 1 + 1. 

q p N 

As (V/9 n ) nS N may be bounded in C([0,T];L p ) and, by embedding, (dtu n ) n £N may be bounded 
in L 1 ([0, T]; L ), in terms of Sobolev norms only, it is thus possible to get a bound of VII™ in 
L 1 ([0, T]; L p ) in terms of the initial data. The rest of the proof goes by the steps that we used 
before. 

5. The proof of Theorem 02 

For T > 0, let us introduce the set Ft of functions (a, u, VII) such that 

a e C b ([0,T] x R N ), Da G C([0, T]; B^ 1 ), 

u G C([0, T];Bp >r ), Vn G C([0, T];L 2 ) n ^([O, T];B s p>r ). 

Uniqueness in Theorem [3] stems from Proposition [HJ Indeed, we see that, as p < 4, any solution 
(p, u, Vil) in Ft satisfies u G C([0, T]; VF 1 ' 4 ) (according to Proposition [2] and to the remark that 
follows) and VII G C([0,T];L 2 ). Therefore, using the velocity equation and Holder's inequality, 
we get 

(63) (dtu-f) = -(«• Vn + aVn) G C([0, T]; L 2 ). 
Note that, as u and Vp are in C([0, T];L 4 ), we have 

(64) d t p = -u-VpeC([0,T];L 2 ). 

Now, consider two solutions (p±,u±,VTl±) and (/92, «2, VII2) in Ft, corresponding to the same 
data. Then ([63]) implies that 6u := U2 — U\ belongs to C 1 ([0, T]; L 2 ) while ([6^]) guarantees that 
Sp := P2 — pi is in C 1 ([0, T]; L 2 ). So Proposition [8] applies and yields uniqueness. 

Let us now tackle the proof of the existence part of the theorem. We claim that if we restrict 
our attention to solutions which are Ft then the assumptions of Lemma [3] are fulfilled so that 
it suffices to solve System ()18j) . Indeed, it is only a matter of checking whether Qu is in 
C([0,T];L 2 ). Applying Q to the velocity equation of CEHD, we get 



a r\.. r\ £ s~\t _ i r\t„. 
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From the assumptions on /, the definition of Ft and the fact that Q maps L 2 in L 2 , we see 
that the first two terms in the right-hand side are in C([0,T]; L 2 ). Concerning the last term, we 
just use the fact that, as pointed out above, u ■ Vu belongs to C([0, T];L 2 ) so and Q(u • Vu), 
too. 

Let us now go to the proof of the existence of a local-in-time solution for (|18p under the 
assumptions of Theorem [3l Compared to Theorem [I] the main change is that we do not expect 
to have u G C([0,T]; L 2 ) any longer (i.e. the energy may be infinite). However, as the pressure 
satisfies 

-div (oVII) = div (u • Wu) - div /, 
Lemma [2] will ensure that VII G C([0, T];L ) anyway if u ■ Wu belongs to C([0,T]; L 2 ) . In 
view of Proposition [21 Remark [2] and Holder's inequality, this latter property is guaranteed by 
the fact that u G C([0, T\;B* )T ) for some p < 4. 

Once this has been noticed, one may use the same approximation scheme as in Theorem[TJ we 
first set (a°,u°,Vn°) := (cto,uo,0). Next, we assume that (a™, u n , VIP) has been constructed 
over R + , belongs to the space Ft for all T > and that there exists a positive time T* such 
that (|19p is satisfied for all t G [0, T*} and, for suitable constants Co and C (one can take 
0) = 2C 2 ), 

(65) ||Va n+1 (i)|| B «-i < 2||Va || B |-i for all t G [0, T*], 

(66) U n (t) < AU (t) + Co^A)(||div/|| Ji i ( B.-x ) + (p*A r||Q/|Li (L2) ), 

(67) a*||VH"|| Ltl(B 3 jr) < c(||div/|| LKBrrl) + (p*A o y j\(U n ) 2 + \\Qf\\ L z) cfr 

with A := a* + IIDooIIb--!, Cr (t) := ||«o||B« r + ll/llrJ(B| ir ) and : = \W n (t)\W p 

Arguing exactly as in the proof of Theorem [H we see that if we define a n+1 as the solution to 

d t a n+1 + u n ■ Va n+1 , a"t = a 

then (PD is satisfied for all time, Va n+1 G C(R + ; J B|~ 1 ) and 

(68) ||Va" +1 (t)|| srrl < e C /o ^W^HVaoll^-i. 



So if we assume that T* has been chosen so that 

(69) C [ T U n {t)dt < log 2 

J o 



then a n+1 satisfies ([65]) . 

Next, we take u n+1 to be the unique solution in C(M. + ;Bp r ) of the transport equation ([2j 
As before, the right-hand side of ([27]) belongs to C(IR + ; ) and one may use ([28]) . So finally, 
the existence of u n+l G C(IR + ; -Bp jr ) is ensured by Proposition [5j and we have 

\K + HmB lr <e C S> n ^{\\uo\\B^ 

+ j\- C Io Un ^ d ^((\\a n+1 \\ L ^ + ||Va" +1 || Brrl )||VH"|| B , r + ||/||^ p ) dr). 

Therefore, if we restrict our attention to those t that are in [0, T*] with T* satisfying (|69j) . and 
use Inequality ([68]) . we see that for all t G [0,T*], 

U n+l {t) < 2U (t) + CA j ||Vn n || B| r oIt with A := a* + ||Va || B |-i- 

w 

So if we assume that Co = 2C 2 and that T* has been chosen so that 
(70) C 2 o*A n f T U n (t)dt < I 
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then taking advantage of Inequality (i67|) , we see that u n+1 satisfies {55) on [0,T*]. 

To finish with, in order to construct the approximate pressure IF 1 , we solve the elliptic 
equation (j33l) for every positive time. Recall that we have div (u n+1 ■ Wu n+l ) G -Bp^ 1 an d that 

||div(u n+1 • VVu n+1 )\\ B s-i < C{U n+1 f. 

Next, given our assumptions on (s,p,r) we have B s p r =— > W 1,4 . Therefore, since "P maps L 4 in 
L 4 , one may write 

\\u n+1 ■ VVu n+1 \\ L 2 < \\u n+1 \\ L 4 \\Vu n+1 \\ L 4, 

< C\\u n+1 \\ L 4u n+1 \\ w i A , 

< C(U n+1 ) 2 . 

Therefore, the second part of Proposition [7] ensures that Vn n+1 is well defined in C(M+;L 2 ) n 
L} oc (R+;B^ r ) and that 

1 1 1 1 t , ) < C 1 1 div / 1 1 (jB 3- x } + (1 -h 1 1 .D^ 1 1 1 (jEf ,- 1 ) ) ^ ^* C^""-^ 1 ) 2 -h 1 1 1 1^ ) ^) - 

Taking advantage of Inequality (|65p at rank n + 1, one can now conclude that VII n+1 satis- 
fies ([5ZD. 

At this stage we have proved that if Inequalities ([65]) . (f66|) and (|57|) hold for (a n ,u n ,VII n ) 
then they also hold for (a ra+1 , -u n+1 , VII n+1 ) provided T* satisfies Inequality (|70p . One may 
easily check that this is indeed the case if we set 

(71) T* := S up{t>0/p*A t(u (t) + p*A \\divf\\ Ll{B s~i ) + (p*A 

for a small enough constant c depending only on s, p and N. 

Once the bounds in Ft* have been established, the last steps of the proof are almost identical 
to those of Theorem [TJ Indeed, introducing 

a n (t,x) := a n (t,x) - a (x) and u n (t, x) := u n (t, x) - u (x) - [ f(T,x)dr 

Jo 

and observing that 6u n := u n+1 — u n = u n+l — u n , one can use exactly the same computations as 
before for bounding 6a n , 5u n and Vdl n . As a consequence (a n ,u n , VIP) is a Cauchy sequence 
in C([0, T*]; L 2 ). Next, the bounds ([65]) . ([66]) and ([6?]) enable us to show that the limit is indeed 
in Ft and satisfies (|18p . The details are left to the reader. 

Let us finally establish the continuation criterion. Note that Lemma U] still applies in the 
context of infinite energy solutions. Hence, repeating the proof of Lemma [5] and using the loga- 
rithmic interpolation inequality of [15] yields the result. This completes the proof of Theorem [3j 

6. The proof of Theorem @] 

In this section, we aim at investigating the well-posedness issue of System ([1]) in Holder spaces 
C s (which coincide with the Besov spaces B^ ^ if s is not an integer), and, more generally, in 
Besov spaces of type B^ r . Of particular interest is the case of the Besov space B^ x which is 
the largest one for which Condition (C) holds. 

The main difficulty is that the previous proofs where based on the elliptic estimate stated in 
Proposition [7| which fails in the limit case p = oo. In this section, we shall see that the case of 
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6.1. The proof of uniqueness. Note that in the case p > 4 the solution provided by The- 
orem |4]need not satisfy VII E L 1 ([0, T]; L 2 ). Hence dtu need not be in L 1 ([0,T];L 2 ) and the 
assumptions of Proposition [8] are not satisfied. 

So, in order to prove uniqueness, we shall prove stability estimates in LP rather than in 
L 2 . These estimates will be also needed in the last step of the proof of the existence part of 
Theorem HI 

Consider two solutions (pi,u\,VIli) and (p 2 ,u 2 ,'VI[ 2 ) of ([I]). Let a\ := 1/pi and a 2 := l/p 2 . 
As usual, denote 5a := 02 — ai, 5u := u 2 — u±, VdTl := V1I2 — VIli and 5f := /2 — f\. First, as 
div u 2 = and 

(9 t & + U2 • V& = — <5u • Vai, 

one may write 

(72) ||<5»(t)||xp < ||&(0)|| LP + rilVoilUooHduHipdr. 

•/ 

Next, as 

dtdu + u 2 ■ V5u = Sf — 5u ■ Vu± — 5aVUi — a 2 V5tt, 

we have 

(73) \\Ht)h P < \\6u(P)\\lp + J (ll^/llii' + 1| Vt*i|Uoo j> + 1| VHi |U« H&HiP + Vfilllip) dT*. 
Finally, we notice that V<91 satisfies the elliptic equation 

(74) div (a 2 V<H) = div 5/ - div (SaVUi) - div (du ■ V«i) - div (u 2 ■ V5u). 
The key point here is that, owing to dWu 2 = div<5u = 0, we have 

div {u 2 • V&i) = div (5u ■ Vu 2 ). 

Hence Equality (J74"|) rewrites 

a*A<H = div 5f + div ((a* - a 2 )VdI) - div (&VII1) - div (5u • V(«i + u 2 )) 
so that the L p elliptic estimate stated in Proposition [6] implies that 

a»||V<H||iP < c(\\5f\\ LP + (a* - a»)||V<H||iP + ||Vni|| L <» H&IIlp + ||V(«i + u 2 )\\ l °°\\6u\\ L p). 

If the quantity a*/a* — 1 is small enough then we thus have, up to a change of C, 

a*||vai|| LP < c(\\5f\\ LP + ||vni|| L oo||&|| iP + ||v(«i + « 2 )IIl»II«*IIh')- 

Plugging this latter inequality in (f73|) then adding up Inequality (|72|) . we get 

||(&,(5u)(t)||zP < \\(5a,5u)(0)\\ LP 

+C f (\W\\lv + (UVuilU- + ||V«2|U- + ||Vai|U« + ||Vni 11^)11(52,^)11^) dr. 

Applying Gronwall's lemma yields the following result which obviously implies the uniqueness 
part of Theorem 0J 

Proposition 9. Let (p\,u\, and (p 2 , u 2 , VH 2 ) be two solutions of JT]) on [0, T]xR w such 

that for some p £ (1, 00), 

• 5a := a 2 — a\ and 5u := u 2 — u\ are in C([0,T]; L p ), 

• V<SI:= Vn 2 - Vn x is in ^([O,^ ;£»»), 

and for some positive real numbers a* and a* such that a* < a*, 

a* < 01, a 2 < a*. 

There exists a constant c depending only on p and on N such that if 

_* _ ^ 
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and if for all t £ [0, T], 

V{t) := /'(llVuillioo + ||V-u 2 |U- + IIVoiIIloo + HVniHioo) dr < oo 

JO 

t/ien i/ie following inequality is satisfied: 

\\(&z,da)(t)\\ LP <e cv ^(\\(&i,5u)(0)\\ LP +C f ' e- cv ^\\8f(r)\\ LP dr) . 



o 

6.2. A priori estimates. Here we assume that (p, u, VII) is a solution to (pQ) on the time 
interval [0, T] with the B^ r regularity. We want to show that if T has been chosen small 
enough then the size of the solution at time t < T is of the same order as the size of the data. 
First, it is clear that we have 

a* < a < a* . 

Moreover, one may write thanks to Proposition [5j 

(75) l|Va(*)|| B *-i < e Cr ^ N ^*[|Vao[| B .-i ) 



and for the velocity, we have, as in the case p < oo, 

\Ht)\\B^ T < e c/o<INb - dT (lKI| B|o>r 

(76) + /^^^"""^^(ll/ll^, + Ikll^JlVnllB^) dr). 
Note that applying standard IP estimates for the transport equation yields 

(77) \Ht)\\ LP < \\u \\lp + [ \\f\\Lp dr + a* [ \\VU\\ LP dr. 

Jo Jo 

As Propositions [6] and [7] fail in the limit case p = oo, in order to bound the pressure, we have to 
resort to other arguments. Now, dividing the velocity equation of (PQ) by p and applying div , 
we get 

(78) aAU = div / - div (u ■ Vw) + div ((a - a) VII) with a:=l/p 
and, by virtue of the Bernstein inequality, we have 

l|vn|| BSoir < HA^vnii^ + iKid-A^vnii^, 
< c||vn|| iP + ||(id- A_i)vn|| BSoiP . 

On the one hand, in order to bound the L p norm of VII, we simply apply the standard LP 
elliptic estimate (see Proposition [6]) to (fT8j) . We get 

a||VII|| j> < c(\\ Q/||iP + || u\\ lp ||V«||loo + ||a — a||L°° || VII | 

Hence, if a* /a — 1 is small enough then 

(79) 5||Vn||iP < c(\\Qf\\ LP + ||«||£p||V«||r« 



LP 



On the other hand, for bounding the high frequency part of the pressure, one can use the fact 
that Operator V(— A) _1 (Id — A_i) is homogeneous of degree —1 away from a ball centered at 
the origin, hence maps -B^" 1 in B^ r (see e.g. [2], Chap. 2). Therefore we have 

a[|(Id-A_i)Vn[| BSoir < Ca||An|| B; -i, 

< C(||div/|| B? -i + ||div(u • Vu)|| B £-j + ||div((a- a)Vn)|| BS -i). 

In order to bound the second term, one may combine the Bony decomposition and the fact that 
divu = 0. This gives 

div (u ■ Vu) = J2( 2T di^ d j ui + diR{u j ,8^)). 
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Thus applying Proposition El we may write 

[|div(it • Vu)\\ B s-i < C\\Vu\\l°°\\u\\b^ iT - 
Finally, as B^ r is a Banach algebra, we have 

||div((s-a)vn)|| BS _i<C7||o-a||B So)P ||vn||Bs 0iP . 

Putting this together with (|79p . one may conclude that there exists a constant c such that if 

(80) \\a-a\\L^(B^ ir ) < ca 
then 

(81) a([|Vn||iP + ||Vn|| Bgor ) < C{\\Qf\\ LP + ||div/|| BS _i + ||«||i>nflj )P ||V«||L»). 



Let us assume that T has been chosen so that 

I r s ~ 1 

-° oo , r 



T 

(82) C / ||Vu||„.-i < log 2 







and that the initial density is such that 

-n c — 

\\ a o - a r S -i < - a. 

Then ([80]) is fulfilled and, combining Inequalities (f76j) , (f77|) and ([8T|) , we get 

tf(t) < 2C7o(*) + Cp||a || B|o>r f(\\Qf\\ L p + Hdiv/llaj-i + U 2 ) dr 

with 

U(t) := \\u(t)\\ LPr]B s and C/ (i) := HuoH^ae^ r + / \\f\\LPr\B« r dT. 

Jo 

It is now easy to find a time T > depending only on the data and such that both Condition 
(1821 and 

U(t) < 4C/o(t) for all f G [0,T] 

are satisfied. 

6.3. The proof of existence. This is mainly a matter of making the above estimates rigorous. 
We have to be a bit careful though since the data which are considered here do not enter in the 
framework of Theorems Q] and [3l 

As a first step, we construct a sequence of smooth solutions. In order to enter in the Sobolev 
spaces framework, one may proceed as follows. 

For the density, one may consider p$ := p+S n (4>(n~ 1 -)(pQ—'p)) where is a smooth compactly 
supported cut-off function with value 1 on the unit ball of R . Obviously, Pq — ~p is in H°° 
and converges weakly to po — ~p when n goes to infinity. In addition, by using the fact that <f> is 
smooth and that B^ r is an algebra, one may establish that there exists some constant C such 
that for all n £ N, 

||Po - Pl|fl&,, P < C1Ipo-p||b£, iP - 

Similarly, for the velocity, one may set Uq '. — S n {^(f){ji ^*)Tio) and for the source term, f n \ — 
ot n *t (SV^fK^ -1 ) • /)) where the convolution is taken with respect to the time variable only and 
(a„) ng N is a sequence of mollifiers on R. 

Applying Theorem Q] thus provides a sequence of continuous-in-time solutions with values in 
H°°, defined on a fixed time interval. Then applying the above a priori estimates, it is easy 
to find a time T independent of n for which the sequence (p n , u n , VIP)„ g N is bounded in the 
desired space. 

For proving convergence, one may take advantage of the stability estimates in L p . The proof 
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6.4. A continuation criterion. This paragraph is dedicated to the proof of the following 
continuation criterion: 

Proposition 10. Assume that s > 1 (or that s > 1 if r = 1). Consider a solution (p, u, VII) 
to ([1]) on [0,T[xM> N such that for some p G (l,oo) we have 

• P GC([0,T);B^ r ), 

• ueC([0,T);B™ r nLP), 

. VII G L 1 ([0,T);B% Oir n]J'). 
There exists a constant c depending only on N and s such that if for some ~p > we have 

sup \\p(t) — p\\b s < cp and / | Vu| £<x> dt < oo 
o<t<T °°' r j0 

then (p, u, VII) may 6e continued beyond T into a B^ Xir solution of ([I]). 

Proof. Applying the last part of Proposition [5] and product estimates to the velocity equation 
of © yields for all t G [0, T), 

||n(t)||^ r <e C /oll v «ll^^f||n ||^ r 



+ jT e" C ^ IIV«D^ (||/|| B ^ + ||a|| BSojr || Vn|| flSoir ) dr) . 



Let us bound the pressure term according to Inequality (i8Tj) . Combining with (1771) and (fT9l) . we 
eventually get 

II II ^ C 1 f II Vnll r oo / 

IK*)l|i>nB' 0ir < e Jo ii ii* ( 

e~ C $a H Vu lli°° dr ^^^^(ll/lliPnB^ + |MI-W"iBg 0>r ||Vu||L°>>) dr\ 

So applying Gronwall's lemma ensures that u belongs to L°°([0, T);L p f\Bl r ). From this point, 
completing the proof is similar as for the previous continuation criteria. □ 

Remark 4. As in the B% framework, an improved continuation criterion involving IIVuIIdo 

instead of \\S/u\\l°° may be proved for the B s oor regularity, if s > 1. The details are left to the 
reader. 

Appendix A. Commutator estimates 

Here we prove two estimates that have been used for estimating the pressure. The first result 
reads: 

Lemma 6. Let (s,p,r) satisfy Condition (C). Let q be in (— 1, s — 1]. There exists a constant 
C depending only on s, p, r, q and N such that for all k G {1, • • ■ , N}, we have 

\\d k [a, A q ]w\\ L p < Cc g 2~ 9? [|Va[| B «-i||i/;||.B* T . for all q > -1 

With ||(Cg)g>_l||^ = 1. 

Proof. We follow the proof of Lemma 8.8 in [TT]. Let a := a — A_ia. Taking advantage of the 
Bony decomposition (|6|), we rewrite the commutator aa^ 

(83) d k ([a, A q ]w) = d k ([T~, A q ]w) + d k T' AqW a - d k A q T' w a + fo[A_ia, Ag\w . 

From the localization properties of the Littlewood-Paley decomposition, we gather that 

R\= Yl d k ([S g >- 1 a,A 9 ]Aq,w). 

\q'-q\<4 
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Note that is spectrally supported in an annulus of size 2 q . Hence, combining Bernstein's 
inequality and Lemma 2.97 in [2], we get 

||-R g ||iP<C 12 \\VSq>-{a\\L°°\\Aq>w\\LP, 
\q'-q\<4 

whence for some sequence (c q ) q >-i in the unit sphere of i r , 

(84) \\R\\\ L y < Cc g 2-^||Va||ioc||^|| B|r . 

To deal with R^, we use the fact that, owing to the localization properties of the Littlewood- 
Paley decomposition, we have 

R l= 12 9 k (Sq'+2A q wA q/ a). 

q'>q-2 

Hence, using the Bernstein and Holder inequalities and the fact that a has no low frequencies, 

ll-Rgll-LP < C ^2 \\Sq'+2AqW\\L°°\\Aq''Va\\LP, 
q>>q-2 

< C2-^2 q ^ +1 - s) J2 2^')( s - 1 )(2^"f ) ||A gU ;|| L o )(2'?'( s - 1 )||A g ,VS|| L p). 

q>>q-2 

Therefore, by virtue of convolution inequalities for series and because N/p + 1 — s < 0, 

(85) \\R 2 q \\Lv < C7cg2 -5? ||Vo|| B .-i||«j|| ,_N. 

Next, Proposition [3] ensures that, under the assumptions of the lemma, the paraproduct and the 
remainder map B^ r x B* in B^ r . As moreover we have 

(86) ||a|| B - r < C\\Va\\ B s-i, 
one may conclude that 

(87) \\R 3 q \\ L p < Cc g 2-^\\Va\\ B s-i\\w\\ Blr . 

Finally, as the last term R^ is spectrally localized in a ball of size 2 q , Bernstein's inequality 
ensures that 

\\R%p < C2 q \\[A_ 1 a,A q ]w\\ LP . 
Then, resorting again to Lemma 2.97 in [2], we get 

(88) \\R 4 q \\ LP < Cc q 2-^\\Va\\ L oo\\w\\ Blr . 

Putting Inequalities ([84]) . ([85]) . ([57]) and ([88|) together and using Proposition [2] completes the 
proof of the lemma. □ 

Lemma 7. Let ^ > and 1 < p, r < oo. There exists a constant C such that 

\\d k [a,A q ]w\\ LP < Cc q 2~ qc; (\\Va\\L°°\\w\\ B s r + \\w\\ L oo ||Va||_B^ r J for all q > -1 

With ||(Cg)g>_l||^r = 1. 

Proof. We use again Decomposition ([83]) . We have already proved in ([83]) and ([88]) that R^ and 
Rq satisfy the desired inequality. Concerning R^, recall that 

\\R 2 q \\LP<C ^2 \\S q '-iA q w\\L°°\\A q tVa\\ L p, 

q'>q~2 

whence 

\\R 2 q \\ L p < C2^ J2 2{q ~ q ' ] " lklk-2 9 ' ? ||A ? /Va|| iP . 

q'>q-2 
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According to Proposition [3j we have 

\\T^a\\ B ,+i < C||w|Uac||5|| B t+i. 

Hence, as ||a|| B< ;+i < C||Va|| B ^ r , the term satisfies the required inequality. □ 

Appendix B. A Bernstein-type inequality 

Lemma 8. Let 1 < p < oo and u G L p such that SuppS C {£ G 1^ / R\ < |£| < -R2} / or some 
real numbers R\ and R2 such that < Ri < Ri- Let a be a bounded measurable function over 
M. N such that a > a* > a.e. There exists a constant c depending only on N and R2/R1, and 
such that 

(89) ca*^^-Jfl? J \u\ p dx < (p - 1) J a\Vu\ 2 \u\ p ~ 2 dx = - J div (aVu) \u\ p ~ 2 u dx. 



Proof. The case a = 1 has been treated in [9] and readily entails the left inequality in (189D for 
one may write, owing to the case a = 1, 



ca*R 2 J \u\ p dx < p 2 J a*|V-u| 2 |u| p 2 dx < p 2 J a\Vu 



\ 2 \u\ p - 2 dx. 

Let us now justify the right equality in (|89p . In the case p > 2, it stems from a straightforward 
integration by parts. 

Let us focus on the case 1 < p < 2 which is more involved. Smoothing out a if needed, one 
may assume with no loss of generality that a is in W 1)0 ° . Let T e (x) = \/x 2 + e 2 for x G R and 
e > 0. We have 

-/ div(aVn)(T £ (u)) p - 1 r e '(u)dx = (p - 1) / a\Vu\ 2 \T'(u)\ 2 (T £ {u)) p ~ 2 dx 

+ [ a\Vu\ 2 T" {u)(T £ {u)) p - 1 dx. 

In view of the monotonous convergence theorem, 

lim / a\Vu\ 2 \T'(u)\ 2 (T £ (u)) p ~ 2 dx = f a\X7u\ 2 \u\ p ~ 2 dx G W + . 
e^oJ RN 1 J r n 1 111 

Next, we notice that 

(90) |div(aVu)(T £ (u)f- 1 T;(u)| < M^ldiv (oV«)|. 

Now, as a G VF 1 ' 00 and u is a smooth function with all derivatives in L p (owing to the spectral 
localization), one may write 

div (aVu) = aAu + Va • Vit, 

hence div (aVu) is in L p and the right-hand side of (|90f) is an integrable function. So finally 
Lebesgue's dominated convergence theorem entails that 

lim / dw(aVu)(T £ (u)) p - 1 T'(u)dx= [ u\u\ p ~ 2 dw (aVu) dx. 
£^0 Jrn J r n 

Therefore 

(91) (P-1)/ a|Vu| 2 |u| p_2 dx < - / uM p - 2 div (aVu) < 00. 

JR N JR N 

In fact, equality does hold. Indeed, whenever 1 / 0, the term T £ (x)T £ (x) p ~~ l tends to when 
e goes to and 

^wr 1 = i*r 2 — (g/x)2 2 , < \ x r 2 . 

(l + (e/x) 2 ) 2 5 

Therefore, as, according to (|9ip . the function |Vii| 2 |u| p_2 is integrable over R , we get 



lim / a|Vn| 2 r"(n)(r e (n)) p - 1 (ix = 0. 
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On the other hand, as u is real analytic, 

/ a|Vu| 2 ^'(«)(T e («)) p - 1 dx = e p - 2 f a\Vu\ 2 = 0. 

Ju=0 Ju=0 

□ 
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